CONSTRAINING TOPOLOGY WITH THE CMB 
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Abstract 

We present a new data analysis method to study rectangular "small universes" 
with one or two of its dimensions significantly smaller than the present horizon (which we 
refer to as T^- and T^-models, respectively). We find that the 4 year COBE/DMK data 
set a lower limit on the smallest cell size for T^- and T^-models of 3000/i~^Mpc, at 95% 
confidence, for a scale invariant power spectrum (n=l). 



1 Introduction 



In the past few years, mainly after the discovery of CMB anisotropics by COBE/DMR pd| , 
the study of the topology of the universe has become an important problem for cosmologists and 
some hypotheses, such as the "small universe" model |Q, have received considerable attention. 
From the theoretical point of view, it is possible to have quantum creation of the universe with 
a multiply-connected topology JlBl. From the observational side, this model has been used to 
explain the "observed" periodicity in the distributions of quasars [0 and galaxies 0. 

Almost all work on "small universes" has been limited to the case where the spatial sec- 
tions form a rectangular basic cell with sides Lx,Ly,Lz and with opposite faces topologically 
connected, a topology known as toroidal. The three-dimensional cubic torus is the simplest 
model among all possible multiply-connected topologies, in which all three sides have the same 
size L = Lx = Ly = Lz- In spite of the fact that cubic T^-model has been ruled out by CO BE 



results [0, H, |T3|, |T3[, the possibility that we live in a universe with a more anisotropic 
topology, such as a rectangular torus T^, is an open option that has not been ruled out yet. 
For instance, if the toroidal model is not a cube, but a rectangle with sides Ly ^ 

and with one or two of its dimensions significantly smaller than the horizon Rh (= 2cHq ), 





Figure 1: Simulated slcy maps for tlie T^- and T^-models and tlieir S'-maps. (A) T^-model witli 
dimensions (L^, Ly, L^) = Rh (3,3,0.3); (B) T^-model with dimensions (Lx, Ly, L^) = i?_ff(0.3,0.3,3); 
(C) S'-map of the T^-model shown in (A); (D) S-map of the T^-model shown in (B). Both models 
are plotted in galactic coordinates and have a scale invariant power spectrum (n=l). 



this small rectangular universe cannot be completely excluded by any of the previous analyses: 
constraints from the DMR data merely require that at least one of the sides of the cell be larger 
than Rh- 

if the rectangular T^-universe has one of the cell sizes smaller 



As pointed out by M and [13 



than the horizon and the other two cell sizes are of the order of or larger than the horizon (for 
instance, <C Rh and L^^Ly > Rh), the values of ST/T are almost independent of the z- 
coordinate, i.e., the large scale CMB pattern shows the existence of a symmetry plane formed 
by the x and y-axes; and if it has two cell sizes smaller than the horizon and the third cell 
size is of the order of or larger than the horizon (for instance, Lx,Ly <^ Rh and Lz ^ Rh), 
the temperatures 6T/T are aproximately independent of both x and y, i.e., the CMB pattern 
shows the existence of a symmetry axis: values of 6T/T are almost constant along rings around 
the z-axis. We call the former case a T^-model because the spatial topology of the universe 
becomes just in the limit Lx,Ly —>■ oo with being fixed. The later case is denoted a 
T^-model for the same reason (the corresponding limit is ^ oo with Lx,Ly being fixed). 
See Figures ITJA (upper left) and |l|B (upper right). 

Our goal is to show that the COBE /J^MK maps have the ability to test and rule out T^- 
and T^-models. We use a new approach to study these models in which we constrain their 
sizes by looking for the symmetries that they would produce in the CMB, obtaining strong 
constraints from the 4 year COBE /'DMR data. 



2 The method 



The analysis of T^- and T^-models is not an easy task, since there are infinitely many 
combinations of different cell sizes and cell orientations. In order to study these models, we 



choose a statistic in which we calculate the function S{ni) defined by 

.... 1 ^- [f(ft.)-f(MF 

where Npi^ is the number of pixels that remain in the map after the Galaxy cut has taken place, 
hij denotes the refiection of hj in the plane whose normal is fij, i.e., 

Uij = hj -2{hi- hj)hi, (2) 

and o"(nj) and a{hij) are the r.m.s. errors associated with the pixels in the directions fij and hij. 
S{hi) is a measure of how much refiection symmetry there is in the mirror plane perpendicular 
to rij. The more perfect the symmetry is, the smaller S{hi) will be. When we calculate S{ni) 
for all 6144 pixels at the positions fij, we obtain a sky map that we refer to as an S'-map. This 
sky map is a useful visualization tool and gives intuitive understanding of how the statistic 
S'(fij) works. 

In order to better understand S'(fij), we first consider the simple model of a T^-universe 
with Lz <^ Rh- For this specific model, the values of 6T/T are almost independent of the 
z-coordinate, so we have almost perfect mirror symmetry about the xy-plane or, in spherical 
coordinates, 6T/T{6,(f)) ^ 5T/T{tt — 6,(f)). When ft, points in the direction of the smallest 
cell size {i.e., in the z-direction), we have S'(fii) ~ 1; otherwise, S'(fii) > 1. An S'-map for a 
T^-model {L^, Ly, L^) = -R//(3,3,0.3) can be seen in Figure |1|C (lower left). Notice in this plot 
that the direction in which the cell is smallest can be easily identified by two "dark spots" at 
hi z and fij ~ — z. For T^-models, the only difference will be that in the place of the two 
"dark spots" , we have a "dark ring" structure in the plane formed by the two small directions. 
See Figure |Ip (lower right), an S'-map of the T^-model {L^, Ly, L^) = -Rh(0.3,0.3,3). 

From the two S'-maps, we can infer two important properties: first, the direction in which the 
S'-map takes its minimum value, denoted So, is the direction in which the universe is small. For 
a large universe such as Lx, Ly, ^ Rh-, the S'o-directions obtained from different realizations 
are randomly distributed in the sky. Secondly, the distribution of S'o-values changes with the 
cell size, i.e., as the universe becomes smaller, the values of So decrease. From the definition of 
the S-map, it is easy to see that the value of So is independent of the cell orientation. In other 
words, if we rotate the cell, we will just be rotating the S-map, leaving its minimum value S'o 
unchanged. 

From here on, we will present our results in terms of the cell sizes Rx, Ry and Rz, usually 
sorted as Rx ^ Ry ^ Rz and defined as Rx = Lx/Rh, Ry = Ly/Ru and Rz = Lz/Rh- We 
remind the reader that the results are identical for all six permutations of Rx, Ry and Rz. 



(1) 



3 Data Analysis 

If the density fiuctuations are adiabatic and the Universe is spatially fiat, the Sachs- Wolfe 
fiuctuations in the CMB are given by |jTO| 

where r is the vector with length Rh = 2cHq^ that points in the direction of observation 
{0,(1)), Hq is the Hubble constant (written here as lOO/i km s~^ Mpc~^) and (5k is the density 
fluctuation in Fourier space, with the sum taken over all wave numbers k. 
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Figure 2: Cumulative probability distribution of 6*0 for T^- and T^-models obtained from Monte 
Carlo simulations. A, left plot: Simulations for T^-universes with dimensions {Kxi Ry, Rz) = (0.5,0.5,3) 
or dot-dashed line, (0.6,0.6,3) or dashed line, and (0.7,0.7,3) or dotted line. B, right plot: Simulations 
for T^-universes with dimensions {Rx, Ry, Rz) = (0.5,3,3) or dot-dashed line, (0.6,3,3) or dashed line, 
and (0.7,3,3) or dotted line. In both pictures the model {Rx, Ry, Rz) = (3,3,3) is represented by a 
solid line, S^^^^ = 2.59 (vertical straight line) and the horizontal solid lines indicate the confidence 
levels of 95%, 90% and 68% (from top to bottom). 



In a Euclidean topology the universe is isotropic, and the sum in (^ is normally replaced 
by an integral. However, in a toroidal universe this is not the case. In this model, only wave 
numbers that are harmonics of the cell size are allowed. As a result, we have a discrete k 
spectrum |P, |15| 

Rh \Rx Ry Rz ) ^ 

where p^,, Py and are integers. 

From equation (|^), we can construct simulated skies by calculating 

^(^,0)cx biCos(27r7)+^2sin(27r7)]a^e-(^^®'=^)'/2^n(^,0), (5) 

Px,Py,Pz 

where gi and g2 are two independent Gaussian random variables with zero mean and unit 
variance, 7 = (^^x + + , a = (^jf-^ + (-^j + (^) oc k"^ and n is the spectral index 
of the scalar perturbations. The term Q-(RH&k^f /2 ^presents the experimental beam function, 
where k± is the length of the k-component perpendicular to the line of sight and G is the width 
of the Gaussian beam given by 6 = FWHM/V8 In 2 ^ 0.43 FWHM, where FWHM is the full 
width of the beam at its half maximum. Finally, we model the noise n{9, cj)) at each pixel i as 
independent Gaussian random variables with mean (nj) =0 and variance (nirij) = (Tij^ij 0. 

We generate our simulated skies as standard DMR maps with 6144 pixels for n=l, with a 
Galaxy cut of 20°, FWHM = 7°, the monopole and dipole removed, add noise and normalize to 
(770 = 34.98/iK (the r.m.s. value at 7° extracted from our DMR map, a 4 year combined 53 plus 
90 GHz map with monopole and dipole removed). Fixing a cell size, we construct a simulated 
sky according to (|^), we smooth this once more by 7° and use the statistic defined in (|I]) to 
obtain an S'-map from which we extract its minimum value So- Repeating this procedure 1000 
times, we obtain the probability distribution of So for that fixed cell size. When we repeat this 
same procedure for different cell sizes, we are able to construct plots as shown in Figure S. 
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Figure 3: Grid of cumulative probability distributions of 5*0 for T^- and T^-models obtained from 
Monte Carlo simulations. The thin-shaded, thick-shaded and grey regions correspond, respectively, to 
the models ruled out at 68%, 90% and 95% confidence. 

In Figure (left plot), we show the cumulative probability distribution of 5*0 obtained from 
the Monte Carlo simulations for the cell sizes {Rx, Ry, Rz) = (0.5,0.5,3), (0.6,0.6,3), (0.7,0.7,3) 
and (3,3,3). The horizontal lines indicate the confidence levels of 95%, 90% and 68% (from top 
to bottom). Comparing these curves with the value S^^^ = 2.59 (represented in the plot by 
the vertical straight line), where 5'^*^^ is the So value extracted from our data set, we conclude 
that T^-models with smallest cell sizes R^, Ry ^0.5 can be ruled out at 95% confidence. As the 
second cell size Ry is increased, the curves shift to the left of the T^-models and we can rule 
out T^-models for R^ ^0.5 at a similar confidence level, see Figure ^ (right plot). In this plot, 
we show the cumulative probability distribution of 5*0 obtained from Monte Carlo simulations 
for the cell sizes {R^, Ry, Rz) = (0.5,3,3), (0.6,3,3), (0.7,3,3) and (3,3,3). 

A more complete picture of the cell size limits is obtained when we construct a two- 
dimensional grid for different values of the cell sizes {Rx, Ry, Rz) with Rz = 3.0 and 0.2 < 
Rx,Ry < 3.0 (see Figure H). The thin-shaded, thick-shaded and grey regions correspond, re- 
spectively, to the models ruled out at 68%, 90% and 95% confidence. Notice in this plot that all 
contours are almost L-shaped, which means that, to a good approximation, the level in which 
a model {Rx, Ry) is ruled out depends only on the smallest cell size, Rmin = min{Rx, Ry}. We 
see that Rmin ^0.5 at 95% confidence. 

4 Conclusions 

We have shown that the COBE/DMR maps have the ability to test and rule out and 
topological models. We have presented a new statistic to study these anisotropic models which 
quantifies the "smallness" of a sky map in a single number. So, which is independent of the 
cell orientation, is precisely sensitive to the type of symmetries that small universes produce, 
is easy to work with and is easy to interpret. 



From the COBE/DMR data, we obtain a lower limit for T^- and T^-models of Rx ^0.5, 
which corresponds to a cell size with smallest dimension of L=3000/i^^Mpc. This limit is 
at 95% confidence and assumes n=l. Since the topology is interesting only if the cell size is 
considerably smaller than the horizon, so that it can (at least in principle) be directly observed, 
these models lose most of their appeal. Since the cubic case has already been ruled out as an 
interesting cosmological model 0, and T^- and T^-models for small cell sizes are ruled out, this 
means that all toroidal models (cubes and rectangles) are ruled out as interesting cosmological 
models. 
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